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Abstract 

The absorption and extinction spectra of sub-wavelength cylinder arrays are shown to present 
two different kind of resonances. Close to the Rayleigh anomalies, the diffr active coupling with the 
lattice periodicity leads to sharp peaks in the extinction spectra with characteristic Fano line shapes 
for both s and p-polarizations. When the material exhibits an absorption line or in the presence 
of localized surface plasmon/polaritons, the system is shown to present resonant absorption with 
wider and symmetric line shapes. For s-polarization our analysis predict a theoretical limit of 
50% of absorption. Interestingly, for p-polarized light and an appropriate choice of parameters, a 
subwavelength cylinder array can present perfect (100%) absorption. 

PACS numbers: 42.25.Bs,42.79.Dj,44.40.+a,73.20.Mf 



1 



The study of the extinction spectra of nanoparticles have drawn much attention recently 
for their potential applications to chemical and biological sensors as well as for Surface 
Enhanced Raman Scattering (SERS) [1-3]. While for a single nanoparticle, the excitation 
of localized surface plasmon/polariton resonances (LSPR) lead to well defined extinction 
line shapes at specific wavelengths [2], for an assemble of nanoparticles the spectra can be 
substantially modified by multiple scattering effects. For absorbing materials, resonant 
radiative diffraction as well as resonant absorption processes could both lead to sharp peaks in 
the extinction spectra. Our main purpose is to understand the relative role of both processes 
in the spectra of nanoparticle arrays. This is an important issue from a fundamental point 
of view but it is also specially relevant for a wide range of applications: By Kirchhoff's 
law [4], the control of absorption by nanostructured materials is equivalent to tailor the 
thermal emission. This is of interest for thermophotovoltaic applications and for the design 
of efficient infrared sources. The recent development of coherent thermal sources [5, 6] has 
stimulated further work on thermal emission from photonic crystals [7-9] or nanoparticle 
arrays [10]. 

In the visible range, plasmon resonances on linear array of metallic particles have been in- 
tensively studied considering the influence of different parameters on the extinction spectra 
[11-13]. Enhanced absorption by dielectric particles resulting of phonon-polariton reso- 
nances has also been demonstrated in the infrared [14]. Peaks in the extinction spectra are 
commonly associated to the excitation of LSPR modified by the diffractive couphng with the 
lattice periodicity [13, 14]. For cylinder arrays, s-polarized electromagnetic radiation (with 
the electric field parallel to the cyhnder axis) can not excite any surface plasmon. However, 
the extinction spectra of subwavelength cylinder arrays of non- absorbing dielectric materials 
present sharp peaks even for s-polarization [15]. 

In this letter, we study the extinction spectra of sub-wavelength cylinder arrays. As we 
will see, the spectra present two different kind of resonances. Close to the Rayleigh anoma- 
lies, the diffractive coupling with the lattice periodicity leads to sharp peaks in the extinction 
spectra with characteristic Fano line shapes for both s and p-polarizations. These geometric 
resonances, associated to radiative coupling in absence of absorption, may lead to resonant 
absorption for appropriate grating parameters. We analytically derive the conditions for 
resonant absorption as a function of the geometry and material's parameters. Another kind 
of absorption resonances, with wider and symmetric line shapes, appears when the material 
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exhibits an absorption line or in the presence of LSPRs. In contrast with geometric reso- 
nances, these absorption peaks are almost isotropic with a weak dependence on the angle of 
incidence and geometry. We will demonstrate that for s-polarization there is a theoretical 
limit of 50% of absorption. Interestingly, we will show that, for p-polarized light and an ap- 
propriate choice of parameters, a subwavelength cylinder array can present perfect (100%) 
absorption. 
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FIG. 1: (s-polarization ) Extinction in s-polarization in a map of frequency (D/X) versus the 
transversal momentum of the incoming radiation (Qq = 27rsin0/A), for an array of SiC nanocylin- 
ders with period D = 4.5/im and radius a = 0.2fj,m. Around D/X =0.36 (A = 12.5/im), there is an 
isotropic extinction peak due to the absorption line of SiC. The inset shows the extinction spec- 
trum, which exhibits a typical Fano line shape, for an incident angle 6 = 15° around the geometric 
resonance (close to the first Rayleigh frequency, i.e. D/X = 1 — QoD/2'k). 

Let us consider an infinite set of parallel cylinders with their axis along the z-axis (see 
inset in Fig. 1), relative dielectric constant e = e' + ie" and radius a much smaller than the 
wavelength. The cylinders are located at r„ = nDu^ = x„Uj. (with n an integer number). For 
simplicity, we will assume incoming plane waves with wave vector ko-Lu^ ( i.e. the fields do 
not depend on the z-coordinate) , with k = cj/c and ko = sin 6 u^+k cos 6 Uy = QoUx+qoUy. 
The reflectance, R and transmittance, T of the cylinder array can be calculated by using 
a standard multiple scattering approach in the dipolar approximation [16], as described in 
Ref. [15]. We define the absorptivity. A, as: A = 1 — R — T and the normalized extinction, 
E, as the ratio between scattered plus absorbed powers and incoming power (notice that, 
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below the onset of the first diffraction beam, the extinction is simply given by the sum of 
absorption and specular reflection, E — A + R). 

Let us first consider the simpler case of s-polarized electromagnetic waves. For a single 
sub-wavelength cylinder, the polarizability is given by 

^ ={C(6'-l) + ...}-z|l + C6" + ---} (1) 



where = 7r(A;a)^|e — Ip. It is worth noticing that (^—Q{l/{k'^azz)} ~ i) so, in 

absence of absorption (e" = 0), the expression above is consistent with the optical theorem 
[15, 17, 18]. Multiple scattering effects, due to the presence of the other scattcrcrs, can be 
included in a renormalized polarizability, S^z [15, 17]. Generalizing the results of Ref. [15] 
to include absorption we find 



(2) 



where Gb is the depolarization term [15, 17, 19], defined as Gb = limr_+ro [G'(r) — Go{r,rQ)], 
being G and Gq the Green function of the periodic array and the free-space respectively, 
and Q'lG} = S{G(0)}. The absorptivity and normalized extinction can now be written in 
terms of azz as: 

Dqo ^ ' 

Eis) = A(^) + ( Q{G} - . (4) 

In order to illustrate the main physics involved in the different reso- 
nant phenomena, we will consider a typical dielectric constant given by: 
e — e^o — uj — i^u) (0;^ — u — i^uj . This is a standard form for a polar material 
(Lorentz model), for a metal, the Drude model is recovered taking cut — and cul — 
Silicon carbide (SiC) nanowires provide a simple model system: its dielectric constant is 
given by this form with the following parameters [20]: 600 = 6.7, ujl — 1.825 10^^ rad.s"^, 
LOT = 1.494 10^^ rad.s"-*^, 7 = 8.9662 lO^-*^ rad.s"-*^. Figure 1 displays the extinction in the 
infrared range, in a map of frequency (-D/A) versus the transversal momentum of the 
incoming radiation {Qo), for an array of SiG cylinders (with D — 4.5 /im and a — 0.2 fim). 
The extinction spectra show two different kinds of resonances: i) close and below the 
Rayleigh frequency cu ^ u}[ ^ — c\Qo — 27r/D\ the spectra present a very sharp peak with a 
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strong dependence on the angle of incidence, and ii) a broad peak close to the absorption 
line of SiC {cu ~ cut) which is almost isotropic. The peaks in the extinction spectra correlate 
with corresponding maxima in the absorption. However, the relative strength of the peaks 
depends on the material and geometrical parameters. This is illustrated in Fig. 2 where we 
have plotted the real part of l/a^z and the imaginary parts of a^z (2a), the absorptivity 
(2b) and the extinction (2c) versus the wavelength for an incident angle 9 — 15° an for 
different lattice constants, D (the inset in Fig. 1 corresponds to a zoom of the bold solid 
line in Fig. 2c). 

The different extinction/absorption resonances resemble the recently discussed "lattice" 
[15, 21] and "site" [22] resonances in absence of absorption. Resonant processes arise when 
the real part of l/{k'^azz) — Gb vanishes. Approaching the threshold of the first propagating 
(diffraction) channel, the real part of Gb goes to infinity as ~ {cuf — a;^)~^/^c/(2D) and can 
compensate exactly the real part of 1/k'^azz giving rise to a "geometric" resonance [15, 17]. 
This can only happen when the real part of the polarizability is positive (e' > 1). As a 
consequence, in s-polarization there is no resonant absorption for metallic cylinders. The 
typical Fano line shape of these resonances is illustrated in the inset of Fig. 1. While it is 
possible to obtain a 100% extinction of the beam ( due to the rcficction resonances which 
appear when the absorption is weak [15] ), the maximum absorption A^ax is limited to a 50 
% (see Fig.2(b,c)). It is easy to show that the highest absorption, Aj^ax — 1/2, takes place 
below the first Rayleigh frequency when 

e" 1 

" 7r(A;a)2|e-l|2 " 2^' ^^-^ 

This is one of the central results of this work. We can also notice that at the Rayleigh 
frequency, cu — cui, the array of cylinders is transparent as the absorption and the extinction 
goes to zero. 

Material or "site" resonances are associated to the zeros of the real part of l/a^^. As 
shown in Fig. 2 (a), it exhibits two zeros, one close to a; ~ a;T the other to a; ~ ujs = 
(eoo^i — ^t) / (coo — !)• The second one corresponds to a small value of ^jja^^}, which leads 
to a weak absorption. The broad resonances in Figure 2(b) and (c) coincide with a maximum 
of Q{azz} and correspond to a resonant absorption of the material. Then, they should not 
be very sensitive to the lattice parameters and order. However, the maximal absorptivity 
value of 1/2 can only be reached at the condition given by eq. 5. Close to the Rayleigh 
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FIG. 2: (s-polarization ) (a): Real part of l/ozz (-) and imaginary part of azz (- -) versus A 
for SiC cylinders. Absorption resonance is due to a zero (— • — ) of ^{1/azz}, which coincides 
with a maximum of Qla^z}- (b): Absorption spectrum versus A for arrays of cylinder of radius 
a = 0.2 nm, for an incident angle 6 = 15° and for several periods D = 4.5 /im (green bold solid 
line), D = 6.67 (red solid line), D = 8.83 /xm (black dotted-dash line), D = 11 fj,m (blue dashed 
line), (c): Extinction spectrum for the same parameters than (b). 

(p-polarization ) (d): Real part of 1/ap (-) and imaginary part of Op (- -) versus A for SiC 
cylinders, (e): Absorption spectrum versus A for arrays of cylinders of radius a = 0.5 fim and 
different periods D = 5.71 /im (green solid line), D = 6.0 /xm (red dashed line), D = 6.15 /xm 
(black bold solid line), D = 6.667 /xm (blue bold dashed line), (f): Extinction spectrum for the 
same parameters than (e). 



frequencies, there is a small blue or red-shifted depending of the value of the real part of Gb 
(see Fig 2a. of ref [15]). 

Let us now consider an incoming wave with the magnetic field parallel to the cylin- 
der axis (p-polarized wave). Following the notation of ref. [15] we can now define 
the renormalized polarizabilities of the effective dipoles pointing along the x and y axis: 



6 



Oixx = (l + ttj^j^^^Gft) a^x, ayy = {l + ayydlG},y ayy, being d^yGb the depolarization 
terms due to the components x and y of the electric field scattered by all the dipoles (except 
the considered one) and axx = ocyy = oip the polarizability of a cylinder in p-polarization, 

' ^%?-l) + ■■■]-^{Uc^' + ■■■]■ (6) 



The general expression for the absorption is: 

A^'^ - -^{Ql\^yy? + Ql\^xx?)Cer (7) 

Hence, in p-polarization, the absorption results from the sum of the contributions of two 
dipoles, one pointing in the x-direction (with an effective polarizability oixx) and one pointing 
in the y-direction (with an effective polarizability (iyy). 

At the threshold of the first propagating order, there is a resonant couphng of electric 
dipoles pointing along the y-axis which leads to the divergence of 3?{9^G';,} ^ — (<^i — 
ij^)'~^/'^uj\/ {2cD) at the Rayleigh frequencies (in contrast '^{dyGh} remains finite). Close to 
the first Rayleigh anomaly (c<j < uoi) and providing that Kjl/ap > 0} [15], the absorption 
contains two parts: a resonant part due to the resonance of cxyy (and similar to the case in s- 
polarization) , and the term due to olxx to which we will refer as a "background contribution" . 
For p-polarization, the resonant part has also a maximal value of 1/2 and reaches this 
maximum when 

Depending on the contribution of the non-resonant part, the absorption can reach a value 
higher than 1/2. 

Site resonances for p-polarization are again associated to zeros of the real part of 1/ axx — 
1/ayy — l/cKp. '^{1/ap} exhibits two zeros (see Fig. 2 (d)): one close to a; ~ los, but 
corresponding to a small value of and to e' = -|-1 and thus a weak absorption (for 

SiC, \s = lO.Oyum), the other close to ~ Up (e' ~ —1) which gives rise to ( phonon- 
polaritons or plasmon-polaritons for metals) LSPR (for SiC, Xp = 10.57/im). In contrast 
with the geometric resonances, the LSPR will lead to the resonance of both dipoles. As each 
dipole can contribute up to 1/2, the absorption may reach 100%. 

We have to notice that unless ^{d^Gf,} — ^{dyGf,}, the two effective dipoles will not 
resonate exactly at the same frequency, thus this value of the absorption maximum will 
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FIG. 3: (p— polarization) : Extinction map for an array of SiC cylinders with a period D = 
6.15 fim, a radius of a = 0.5 /um. The localized surface phonon-polaritons resonance, at the 
wavelength A ~ Ap, is isotropic. 

be very sensible to the lattice parameters. Together with this condition , it can be easily 
shown that the other condition to reach the maximal absorption is go = Qo(^ = 45°). We 
represent in Fig. 2(e) and (f), respectively the absorption and extinction spectra of an array 
of SiC cylinders with a radius of a = 0.5 fim and for an incident angle of ^ = 45° and 
several periods. Similarly to s-polarization, there are two kinds of resonances, (1) geometric 
resonance characterized by very sharp peaks at wavelengths close to Rayleigh frequencies, 
(2) a broader double peak at a wavelength A ^ Ap. This double peak corresponds to the 
resonances of the y— and x— dipoles, (red or blue) shifted depending on the sign (+ or -) of 
d^Gb and dyGb- As discussed above, ~ 100% absorption takes place only when dlGb ~ d'^Gb- 
It is worth noticing that in the wavelength range where 3?{l/ap < 0} (i.e. |e| < 1), there are 
no geometric resonances. This explains the anomalous shape of the LSPR extinction peak 
for D = 6.15 fim (black bold solid line in Fig. 2(f)), as the Rayleigh frequency (A = 10.5 fim) 
lies in a range where ^{l/aa} < 0. At this wavelength, we can just observe the (y— dipole) 
Rayleigh transparency dip in the extinction spectrum (see also Fig. 3). These results provide 
a simple and analytical explanation of the anomalous shape of the extinction peaks observed 
in recent numerical simulations [13]. 
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In conclusion, this paper gives a simple analytical method to derive the optical proper- 
ties of sub-wavelength cylinders arrays. We have shown that in s-polarization, absorption 
resonances can absorb up to half of the incident power. The absorption and extinction 
spectra for p-polarization presents a more complex structure due to the contribution of two 
orthogonal dipoles. It is remarkable that, as we have shown, by an appropriate choice of 
parameters a subwavelength cylinder array can become a perfect absorber. We believe that 
our analysis paves a new way for the nanoengineering of chemical and biological sensors and 
photo-thermal devices based on nanoparticle arrays. 
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